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Abstract 

In this work is presented the software OGCOSMO. This program 
was written using high level design methodology (HLDM), that is 
based on the use of very high level (VHL) programing language as 
main, and the use of the intermediate level (IL) language only for the 
critical processing time. The languages used are PYTHON (VHL) 
and FORTRAN (IL). The core of OGCOSMO is a package called 
OGCJib. This package contains a group of modules for the study of 
cosmological and astrophysical processes, such as: comoving distance, 
relation between redshift and time, cosmic star formation rate, number 
density of dark matter haloes and mass function of supermassive black 
holes (SMBHs). The software is under development and some new 
features will be implemented for the research of stochastic background 
of gravitational waves (GWs) generated by: stellar collapse to form 
black holes, binary systems of SMBHs. Even more, we show that the 
use of HLDM with PYTHON and FORTRAN is a powerful tool for 
producing astrophysical softwares. 

Keywords:Computational Physics, Cosmology, Gravitational Waves, Black 
Hole, Structure Formation. 

1. INTRODUCTION 



Around 300 thousand years after the Big Bang, matter and radiation de- 
coupled and the photons could freely travel through the space. Henceforth, 
today we can obtain information about this period from the Cosmic Mi- 
crowave Background Radiation (CMBR). Moreover, about 840 million years 
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after the Big Bang, the Universe was reionized by the first stars. However, 
we have no direct information from the period between 300 thousand years 
and 840 million years. This period is so-called dark age (DA). On the other 
hand, the theory of general relativity predicts the existence of gravitational 
waves (GWs) as perturbations of space-time which propagate at the speed 
of light. The detection of GWs will open a new astronomical window for 
observing the Universe. In particular, they will allow us a deeper understand 
about the DA. 

In this context, we are interested in to describe, by analytical and semi- 
analytical models, some astrophysical and cosmological processes as, for ex- 
ample, the cosmic star formation rate, and its connection with the stochastic 
background of GWs [4J, this kind of study will shed light on the knowledge 
of the process that took place at the end of the DA. 

As such, this modeling framework will provide us with an important aux- 
iliary tool for the reconstruction of the cosmic history. Thus, we decided to 
organize the computational development as a software, it called OGCOSMO, 
whose main characteristics are discussed in this work. 

In section 2 it will be showed a short description of the astrophysical 
and cosmological model adopted. In section 3 it will be presented the main 
characteristics of the software and the type of programing used. In section 4 
we present some results. Finally, in section 5 are the final considerations. 

2. The Model 

In this work we are considering the general theory of relativity with cold 
dark matter and cosmological constant (ACDM model) for background cos- 
mology. 

Here is assumed the hierarchical scenario of structure formation, having 
as base the Press-Schechter (PS) like formalism [6]. The basic idea behind 
this scenario is that the formation of objects like galaxies and galaxy clusters 
occur in the following way: First, when the mean density of dark matter 
perturbation, within a given volume, is larger than a threshold level, 8 C , 
the matter leaves the linear regime and collapses to form small-halo objects. 
These halos become gravitational wells that attract the baryonic matter, that 
is the ordinary matter that form stars and planets. Durant all this process, 
greater halos are formed by fusion of the minors and more baryons fall into 
structures. The star formation starts and black holes grow up by accretion 
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of matter. The complete details, about the all considerations and the main 
results obtained with this model, can be seen in [U [5]. 

3. The Software 

The OGCOSMO software was written using object-oriented programing 
paradigm, that permits to construct codes that are really reusable and clean. 
The main programming language used was PYTHON [7j, that is a very-high 
level language programing, and FORTRAN only for the critical time parts 
of the code. This form of writing codes is called of high level design [3j [2] . 
The external modules used were basically: Tkinter [5] for construction of the 
Graphical Using Interface (GUI), scipy and numpy for numerical methods 

®- 

The core of the OGCOSMO is the OGC Jib, that is a package containing 
the following principal modules, up to now: 

• OGC_cosmo: This module contains the class "CosmoQ" that is com- 
posed by the cosmological background methods, and callbacks, that 
are, comoving distance and comoving volume, relation between time 
and redshift (z), matter density evolution (dark and baryonic mat- 
ters), grow function of matter perturbation, variance of matter linear 
density field, linear extrapolation of the critical density for collapse of 
structures in a given z , comoving volume of dark matter halo. 

• OGC_PS: This module contains the class "PressSchechter(Cosmo)" , 
that inherit the methods of Cosmo class. This class contains the base 
of Press-Schechter like formalism, that is, the functions used for study 
of structure formation as, mass function and numerical density of dark 
haloes, fraction and accretion rate of baryonic matter within structures 
and the cosmic star formation rate - CSFR. 

• OGC_SMBH: In this module is the "SuperBH(PressSchechter)" class, 
that is a class under development and it will be used for the study of 
the evolution of supermassive black holes through of its mass function 
and by the coalescence of these objects [5] . 

4. Results 
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4.1 The OGCOSMO 



In the figure [T] we show the GUI of OGCOSMO. In A are the spaces for 
entry parameters such as the cosmological and those associated with the star 
formation rate. In B are the buttons that start the calculus of the CSFR, 
supermassive black holes mass function and some buttons for creation of 
specific plots as the CSFR (see figure [2]) . 

In the figure [2] is represented the CSFR behavior, whose result can be 
seen by pressing the button "Grafico CSFR" . The Salpeter exponent of initial 
mass function is x = 1.35, the characteristic scale for star formation is r = 
2.5 x 10 9 years and the exponent of the star formation rate n = 1. For more 
details see jlj. 
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Figure 1: GUI of OGCOSMO. In A are showed the entry parameters. In B 
are showed the buttons for running the models and generate some graphics. 




Figure 2: Plot of the cosmic star formation rate behavior. In this case 
was considered fl m = 0.24, Qf, = 0.04, = 0.76 and h = 0.73, for the 
cosmo logical parameters, x = 1.35 and r = 2.5 x 10 9 years. 
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4.2 OGCJib Usage 



In Python a package is basically a hierarchical file directory structure 
which defines a environment that consists of modules. This directory contains 

a file named _Jni py that identifies the directory as a package [lj. In the 

context of this work, this means that if a user does not want to use the GUI 
of OGCOSMO, he (or she) can use only determined methods by importing 
its specific class from a module within OGCJib. For example, an user want 
to call the method age from the class Cosmo of the module OGC_cosmo: 

[1] >>> from OGCJib. OGC_cosmo import Cosmo 

[2] »> MyUniverse = Cosmo(0.04,0.24,0. 73,0.76,6.0,20.0,'. /trabalho') 

[3] >>> Age = MyUniverse. age (5) 

[4] »> print' Age = %3.9e' %Age 

[5] Age = 1.189273236e+09 

In line [1] was imported, from OGC_cosmo module within OGCJib pack- 
age, the class Cosmo. In line [2] was created a new object (MyUniverse) 
and in line [3] was used the method age. The example above show that for 
Q m = 0.24, VL b = 0.04, Ha = 0.76 and h = 0.73, the age of the Universe at 
z = 5isl.2xl0 9 years. The last argument of the class Cosmo ('./trabalho') 
is the directory where the data will be saved. 

This result shows that OGCJib can be used as a framework for fast an 
easy development of astrophysical and cosmological applications for the study 
of the cosmic history, within hierarchical scenario for structure formation. 

5. Final Considerations 

In this work are presented both the initial stage of the OGCOSMO soft- 
ware and the OGCJib package, that have been developed to be an auxiliary 
tool for reconstruction of the cosmic history. Moreover, it is possible for 
an user to have access for a specific method of the OGCJib. This is an 
important feature to show that this package can be used as a framework 
for construction of others softwares. An another aspect is that the program 
can be used as didactic tools in lectures of theory of general relativity and 
cosmology. 
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The next steps will be to write classes to obtain the mass function of 
SMBHs; research of stochastic background of gravitational waves (GW) gen- 
erated by: collapse of stars to black holes, binary systems of SMBH; calculate 
the signal/noise rate for GW detectors (LISA, LIGO, Decigo, BBO). 

Even more, it was showed that the use of PYTHON with FORTRAN for 
high-level design program is a powerful tool for development of astrophysical 
softwares. 
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